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Abstract: We present the laboratory verification of a method for re¬ 

moving the effects of frozen-flow atmospheric turbulence using a Linear 
Quadratic Gaussian (LQG) controller, also known as a Kalman Filter. 
This method, which we term “Predictive Fourier Control,” can identify 
correlated atmospheric motions due to layers of frozen flow turbulence, and 
can predictively remove the effects of these correlated motions in real-time. 
Our laboratory verification suggests a factor of 3 improvement in the RMS 
residual wavefront error and a 10% improvement in measured Strehl of the 
system. We found that the RMS residual wavefront error was suppressed 
from 35.0 nm to 11.2 nm due to the use of Predictive Fourier Control, and 
that the far held Strehl improved from 0.479 to 0.520. 
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1. Introduction 

We present progress towards the implementation of a Predictive Fourier Control (PFC) III] El 
algorithm on sky to correct for frozen flow turbulence in adaptive optics systems. Recently, 
linear-quadratic Gaussian (LQG) controllers Ea have been used for the stable control of 
adaptive optics systems a , and for the suppression of specific frequency vibrations in operating 
adaptive optics systems 00. Predictive Fourier Control is an LQG method which explicitly 
corrects for frozen-flow turbulence crossing the telescope aperture. It diagonalizes that problem 
by working in Fourier space, making the method computationally feasible for high-order AO 
systems. 

Several groups have formulated LQG controllers for the optimal control of adaptive optics 
systems. Gavel and Wiberg a initially developed a Kalman filter for optimal AO control. Le 
Roux et al. a formulated an LQG controller for adaptive optics. Petit et al. a demonstrated 
the use of LQG in the laboratory to suppress vibrations. Sivo et al. a have applied this on sky 
with the CANARY demonstrator. LQG control has progressed from the laboratory to the latest 
generation of astronomical instruments. The recently commissioned Gemini Planet Imager 0 
and SPHERE 0 instruments have both demonstrated the use of optimal LQG control 0 fTOll . 
GPFs AO system uses an LQG controller to suppress vibrations in both tip-tilt and focus 0| . 
Similarly, SPHERE uses an LQG controller to suppress vibrations in only the tip-tilt modes 
Qo). These controllers are implemented with the same LQG framework that we will apply to 



correct frozen flow turbulence across spatial frequencies beyond the tip-tilt regime. 

Our experimental setup provides a tightly integrated way to test the use of an LQG controller 
both in the lab and on-sky. The ShaneAO system with the ShARCs Camera on the Shane 3-m 
telescope at Lick Observatory mmH provides an on-sky instrument. At the Lab for Adap¬ 
tive Optics we have configured an adaptive optics test bench to mimic the hardware and 
physical conditions for ShaneAO. These closely coupled systems allow us to test and develop 
techniques in a controlled lab before proceeding to the instrument on-sky for functional tests. 

This paper describes a laboratory implementation of the full Predictive Fourier Control (PFC) 
algorithm. In Sectionj^we provide a brief overview of the test bench, a full, physical simulator 
designed to test adaptive optics hardware, software and techniques. In Sectionj^we describe the 
PFC algorithm, and the way this algorithm interacts with an operating adaptive optics system. 
In Section|^we describe the implementation of the PFC algorithm in the lab, and the measured 
performance improvement. Section [^discusses the significance of this laboratory demonstra¬ 
tion. 

Our framework, described in Section]^ and diagramed in Figure]^ consists of three compo¬ 
nents: 


1. Phase reconstruction as performed by a standard adaptive optics system. 

2. The Predictive Fourier Controller, which uses a Kalman Filter to provide a best estimate 
of the correction in the future. 

3. Fourier Wind Identification, which provides the state space model to the Predictive 
Fourier Controller. 


Fourier Wind Identification is described theoretically in SectionjTT] and demonstrated in the 
lab in SectionBAl We discuss the theoretical framework of the Predictive Fourier Controller in 


Section 3.2 and demonstrate its operation in Sections 4.2 andpT 


2. Experimental Setup 

Our experiments are performed on a modified and optimized version of the UCSC Laboratory 
for Adaptive Optics (LAO) MCAO/MOAO testbed lfT3ll . In the past, its utility and versatility 
has been shown by simulating MCAO and open-loop MOAO on a 30-m equivalent aperture 
telescope lEl, demonstrating the positive effect of linearity calibrations to Shack-Hartmann 
wavefront sensors during open-loop MOAO operation on a 10-m telescope UB and for devel¬ 
oping wavefront reconstruction and control algorithms for MCAO on an 8-m class telescope 
Da. Currently, we have reconfigured the testbed to match as closely as possible the new adap¬ 
tive optics system on the Shane 3-m telescope at Lick Observatory (ShaneAO) im so that 
the transition from lab to a real system is as seamless as possible. The optical layout of the 
reconfigured system is shown in Figure [T] 

For our purposes, the testbed is operated in single conjugate (SCAO) mode with two de¬ 
formable mirrors (DMs): a low (spatial and temporal) frequency, high-stroke Alpao DM52 
“woofer” and a high-frequency, lower-stroke, Hamamatsu SLM-X8267 spatial light modula¬ 
tor (SLM) “tweeter”. As the test bench does not operate in real time, we treat both mirrors as 
having instantaneous temporal properties and ignore the slower temporal performance of the 
“woofer” mirror. 

Both DMs are conjugate to the equivalent of the atmosphere’s ground layer represented by a 
phase plate mounted on a linear stepper motor. Each phase plate is 160 subapertures across with 
a minimum wind speed of 0.01 subapertures per timestep. The linear stepper can be reliably 
repositioned so that repeated runs can be performed across the same section of a phase plate 




Spatial light 



Fig. 1. Testbed layout for wind prediction tests. The light source is a super-luminescent 
diode laser. The spatial light modulator is a high-frequency deformable mirror. In single 
conjugate mode only one of the Shack-Hartmann wavefront sensors is used. 


and a higher signal-to-noise ratio is achieved by co-adding runs. The size of the phase plate 
ensures that each test run contains 4.4 pupil crossings. 

The woofer’s actuator layout is 4 x 4 across the pupil. The edges of the woofer are pinned, so 
outer actuators have reduced stroke and are not used. This is nearly identical to the conhguration 
of the woofer used in ShaneAO. The SLM has a resolution of 768 x 768, however pixels are 
binned to emulate a 36 x 36 actuator conhguration similar to that of the MEMS tweeter on 
ShaneAO. 

Prior to every set of extended runs (typically once a day), the woofer and tweeter are re¬ 
registered to the wavefront sensor and a new woofer inhuence matrix is generated. Woofer 
edge effects are minimized by performing a ht to individually measured inhuence functions. 
A singular value decomposition (SVD) reconstructor zeroes out potentially problematic (high 
frequency) woofer modes. Each DM has a separately set integrator leak gain (up to 0.99 for the 
woofer and 0.995 for the tweeter), similar to ShaneAO, to allow each mirror to handle its own 
actuator clipping and windup. Detailed telemetry (residuals, DM commands) and conhguration 
is saved during each run to ensure reproducibility. 

The test bench control law mimics the integrators and reconstructors used in ShaneAO, how¬ 
ever the system delay can be arbitrarily shortened or lengthened to match the operating control 
delay measured with ShaneAO or to test the effects of longer or shorter time delays. 

2.7. Modifications to the Experimental Setup to Support Predictive Fourier Control 

In order to use the test bench as a scale representation of the ShaneAO system, and to integrate 
the LQG controller into the test bench software architecture, several changes to the test bench 
were required. 

We hrst re-mapped the pupil of the test bench to the WES lenslet arrays to simulate subape- 
rutres which are scaled equivalents to the 10 cm subapertures found in the ShaneAO system. 
[^The remapping resulted in 36 active subapertures across the pupil, slightly more than are 
found on the ShaneAO system, which will have 30 subapertures across the pupil in its highest 

’ Throughout this paper, we use test bench units, timesteps and subapertures, as our test bench does not operate in 
real time. All of these numbers can be easily scaled to match the 3 m ShaneAO system, an 8 m or a 30 m telescope. 
For the ShaneAO system, we convert from the test bench to an operating frequency of 1 kHz and a subaperture size of 
10 cm. 











































performance configuration. 

As well, we implemented an SVD reconstructor for the system’s high stroke woofer. This 
allows us to eliminate problematic high frequency woofer modes, but to use the natural ac¬ 
tuator basis set for the woofer directly. This also allows us to easily re-measure the woofer 
influence functions for each actuator, recalibrating the system for minor changes in alignment 
and optimizing the system’s long term performance. 

Although the SVD reconstructor for the woofer improved its performance and registration, 
the SLM registration algorithm had to be improved to fix actuator misregistration errors that 
would cause very high spatial frequency content to build up on the mirror. We solved this 
with an improved actuator registration method which uses actuators across the full SLM mirror 
working surface, and automatically hts the pixel locations instead of relying on user input to 
identify actuator pokes. 

The original fiber laser was replaced with a super luminescent diode (SLD) to fix non- 
uniform pupil intensity effects that manifested as an 8 Hz power peak due to laser power 
fluctuations. Over long runs (10,000 iterations or greater), we observe a correlation between 
RMS error and measured Strehl indicating low non-common path errors. The far-held image is 
sharpened using an algorithm that iteratively maximizes Strehl by successively placing Zernike 
modes on the woofer and tweeter. 

We also improved the control loop architecture of the test bench to facilitate the use of the 
Kalman reconstructor. The new architecture allows us to transition quickly between reconstruc¬ 
tors, to perform calibration with one reconstructor (usually a simple Fourier-transform recon¬ 
structor) and then transition to a different reconstructor (the full LQG controller) for operation. 
These modihcations also permit the system to enable, disable and reset integrator states during 
normal operation. 

2.2. Error Budget for the Experimental Setup 

We have developed a consistent error budget for the test bench in its current conhguration, 
based on the error budget presented in Ammons et al. (2010) M- Similar to Ammons et 
ah, we derive many of our error budget terms from the total atmospheric turbulence as meas¬ 
ured by the wavefront sensor’s RMS total wavefront error measured with the phase plate 
held still (Oatmosphere = 180.0 nm). From this measurement, we can derive the htting error 
{Gfi, = Q.Q^lOatmosphere, See QU Section 4.3.2) and aliasing error {OaUas = 0.4oy,>, see CSl 
section 4.3.2 and references therein). The measurement error from the wavefront sensor was 
measured directly on the test bench, with no atmosphere present, and found to be 10.0 nm. 

We explicitly separate bandwidth and delay temporal errors. As the test bench operates in 
discretized timesteps, there is no bandwidth error. The error due to the artihcial delay in the 
system was computed from equation 9.56 in Hardy (1998) ifTTll . 

(^dday=28.4(/,T,)^/' (1) 

where fg is the Greenwood frequency, which can be computed from the wind velocity, and Tj 
is the delay between measurement and correction, which for the experiments described in this 
paper was set to = 2 timesteps. 

To measure the non-common path errors in our Strehl measurement, we applied the extended 
Marechal approximation (see Hardy (1998) equation 4.40 ifTTl . and ifTSll i to the far-held Strehl 
after image sharpening but with no aberrations in the system. After image sharpening, the far- 
held Strehl was 0.80, contributing G^cp = 49.5 nm of error to the system. Examining the error 
budget during experiments suggests that our far-held correction did not remain stable, and 
that an additional 44.0 nm of calibration error exist in the system. All of these error terms are 
presented in Table[2 The RMS far-held wavefront error in Table [^includes terms which are not 


seen by the wavefront sensor. When considering the RMS residual wavefront error as measured 
by the wavefront sensor, we ignore the Mirror fitting, spatial aliasing, static uncorrectable and 
calibration error terms. When measuring the far-held strehl, these terms are included. 


Table 1. Error Budget for the test bench. 


Error Term 

Value 

Wavefront Sensing 

10.0 nm 

Mirror Fitting 

8.8 nm 

Spatial Aliasing 

3.5 nm 

Time Delay 

40.5 nm 

Static Uncorrectable 

49.5 nm 

Calibration 

44.0 nm 

RMS Far-Field Wavefront Error 

78.8 nm 

Predicted Strehl Ratio 

0.471 


For the remainder of this paper, when we discuss RMS residual wavefront error, we are re¬ 
ferring to only those terms which are seen by the wavefront sensor, or to measurements taken 
directly from the wavefront sensor. In contrast, the RMS total wavefront error is the wavefront 
error measured through our simulated atmosphere but with no correction applied. When con¬ 
sidering the wavefront error in the far held, we will always convert to the Strehl ratio using the 
extended Marechal approximation. 

3. Theoretical Framework 

We build on a Predictive Fourier Control (PFC) framework for Adaptive Optics Systems HU El 
that aims to minimize temporal errors, including servo-lag, which manifest as a mis-estimation 
of the current atmospheric phase. This predictive controller uses a slow 0.1 Hz) loop to 
identify frozen How turbulence seen by the AO system, and an LQG controller in the Fourier 
domain to suppress the specihc errors which arise from frozen How turbulence across the tele¬ 
scope aperture. 

The PFC algorithm uses telemetry from the adaptive optics system to identify frozen How. 
Telemetry is recorded over the course of 1-10 seconds and then post-processed. The post¬ 
processing uses the Fourier Wind IdentiHcation scheme (Section[3T] and ifTOl l. as well as a vi¬ 
bration and aliasing analysis to identify specihc temporal frequencies which should be masked 
and eliminated. This information is used to construct a state space model that captures the dy¬ 
namics of the atmosphere. The relationship between the Predictive Fourier Controller and the 
Fourier Wind IdentiHcation is shown in Figure]^ 

This section describes the framework for Predictive Fourier Control, Hrst describing Fourier 
Wind IdentiHcation (Section |3T), then describing the Predictive Fourier Controller (Sec- 
tion |3.2[ ). The application of this framework to the test bench is discussed in Section]^ 

3.1. Fourier Wind Identification 

Using pseudo-open loop phase telemetry from the system, the Fourier Wind IdentiHcation 
(FWI) technique HI identihes atmospheric frozen flow for correction. The Fourier basis set 
provides a convenient way to examine translating frozen-flow, as the individual Fourier modes 
each oscillate at a speciflc temporal frequency. 

The temporal power spectral density (PSD) of each Fourier mode is estimated from the data. 
Individual Fourier modes are split into segments of length S, which are windowed to emphasize 
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Fig. 2. A block diagram showing the dual-loop arrangement used with Predictive Fourier 
Control. The usual AO loop runs at kFlz frequencies, while a separate analysis loop mea¬ 
sures the frozen flow velocities and generates a new Kalman filter for the system. Frozen 
flow turbulence is generally stable for more than 10s of seconds GD, so updating the filter 
on these timescales should not affect the quality of the predictive control. Flowever, even 
if the wind velocity does change before the filter is updated, we have shown (Section [4.3| l 
that the correction does not degrade significantly. 


the middle of the segment, and which overlap with neighboring segments. The overall frame 
rate, fs sets the maximum estimated temporal frequency at /s/2. The length of the segments sets 
the frequency sampling spacing at fs/S. The half overlapped segments increase the signal-to- 
noise in the resulting PSD (for a more detailed discussion of this method, see d). By increasing 
or decreasing S, we can control the signal-to-noise and frequency resolution of the PSDs. 

Once we have created a PSD for each Fourier mode, we can look for temporal peaks which 
are indicative of frozen flow. In Fourier space, temporal peaks will appear with frequencies 
given by 

ft= fxVx + fyVy (2) 

We model each peak as an oscillation at the temporal frequency /, with an added white 
noise broadening term ||2|. Three sample PSDs from the test bench, which show peaks that obey 
Equation (|^, are shown in Figure Each peak corresponds to a potential match to Equation. 
The peaks that appear close to 0 timesteps^* are eliminated as they correspond to the slowly 
varying steady-state errors found in every system. Peaks are fit for each spatial Eourier mode 
separately. 

Using the identified peaks from all the Eourier modes, the EWI algorithm works backwards 
through Equation (|^. The frequencies in Equation Q, when shown on an fx,fy grid, appear 
as a plane in frequency space, with the fx = 0,fy — 0 (piston) term always at 0 timesteps^'. 
Eigure|^ shows the process of matching found peaks in a PSD to a theoretical plane in Eourier 
space. EWI then produces a metric in velocity space that shows the percentage of matched 
peaks at each velocity. Areas with high metric scores are velocities at which frozen flow has 
been detected. 

Once peaks have been identified, they can be formulated into a state-space model. Solving 
the Discrete Algebraic Riccati Equation (DARE) generates a controller that suppresses those 
specific temporal frequencies in its error transfer functions. The gains, leaks and high-pass 
coefficients of the LQG controller are then updated once the frozen flow layers and vibrations 
have been identified. 
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Fig. 3. A temporal power spectral density (PSD) for a three separate Fourier modes 
(^ = 3,6,9 / = 0) which shows clear evidence of frozen flow turbulence. The system 
ran using an LQG controller which did not model any atmospheric frozen flow. Peaks at 
~—0.021,—0.042, and — 0.062 timesteps^* are due to the effect of frozen flow turbulence 
caused by a phase screen mounted on a stepper motor passing in front of the entrance pupil 
to the test bench. The non-zero peak positions are set by Equation 0. The peaks at very 
high frequencies are spatial aliases of the frozen-flow wind. 


MCAOBBnch-2014Q820_222416 



0.0 0.2 0.4 

fx (subapertures”') 



MCAOBench-20140820_222416 



0.0 0.2 0.4 

fx (subapertures”') 


Match 


Possible 


Not Possible 


Fig. 4. Identification of frozen flow turbulence on the test bench using the Fourier Wind 
Identification algorithm. This match shows a frozen flow layer that is moving in the x- 
direction (vx = 0.25 and v,, = 0.0 subapertures/timestep). The left panel shows the detected 
peaks in telemetry. The colorbar corresponds to the temporal frequency of the detected 
peaks. Modes where no matching peak was found are shown as white. The right panel 
shows which peaks match the detected frozen flow layer. It is red for modes which were 
not detected in the telemetry, and white for modes where the detected peak would be too 
close to 0 timesteps^^ to be detected. For this set of telemetry data, 97% of modes where 
peaks could be detected showed evidence of frozen flow. 





























3.2. LQG-Based Predictive Fourier Control 

The wind parameters derived in the Fourier Wind Identification step are then used to calcu¬ 
late an optimal LQG controller for the AO system. The controller uses a Kalman state-space 
framework with an independent state for each frozen flow layer. A derivation of the appro¬ 
priate Kalman filter for predictive control in a variable delay adaptive optics system is done 
analytically and numerically in m and 0, and only a brief overview is provided below. 

Frozen flow in the Fourier domain results in each Fourier mode shifting with a frequency 
given by Equation 0. To apply this shift over a fixed delay between the phase measurement 
on the WFS and the application of the phase to the DM, we use a complex valued gain (ga in 
the derivations from Gl). The complex gain serves to translate each Fourier mode to account 
for the delay between WFS sensing and command application, usually somewhere between 1 
and 2 timesteps. A purely real a corresponds to no frozen flow and a complex a determines 
the shift applied to each Fourier mode, with the complex phas^of a given by 

T 

Phase(a) =-27r(kv;c-|-/vj.)— (3) 

where k and / are the respective Fourier mode numbers, [v;c,V 3 ,] (m/s) is the wind vector for a 
single frozen flow layer, T is the delay between the WFS measurement and the DM command 
application (s), N is the phase grid size {N — 36 for the test bench) and d is the subaperture 
diameter (m). The magnitude of a, |a|, sets the “memory” for an individual layer, and is less 
than 1. Along with a, the noise covariance matrix of the frozen flow layers (including the static 
term) is used to solve the DARE, and yield an optimal controller with ga tuned for each Eourier 
mode. 

The LQG controller can be applied directly with the Kalman Eilter matrix equations. Eor our 
model, it can be simplified into the controller C{z) 

which provides insight into how the prediction works. In Equation Q Q = p -|- cj/ is the total 
layer power, and O’/ is the noise variance of the particular Eourier mode. 

This filter produces transfer functions which show a “notch” at the particular layer frequency, 
as well as good rejection around 0 Hz, the DC term. The theoretical stability of these filters is 
examined in detail in fl] and 0 . 

4. Experimental Results 

4.1. Modeling Transfer Functions 

With the improvements described in Section [ZT] we first demonstrate that the test bench can 
reproduce theoretical transfer functions across a wide range of spatial frequencies for our sys¬ 
tem. To measure the error transfer functions (ETEs), we recorded the wavefront sensor residual 
error measurements for 4096 iterations with the control loop active and closed, and with the 
control loop uncontrolled and open. Eor each telemetry data set, we converted to the Eourier 
basis, and then created periodograms using the technique described in SectionjTT] To compute 
the error transfer function, we took the ratio of the active, closed loop PSD for a given mode 
with the uncontrolled, loop open PSD. An ETE for a single Eourier mode (k = 5,1 = 5) is shown 
in Eigure|^ along with a model for that ETE. 

^This should not be confused with the “phase” of the wavefront. 
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Fig. 5. An empirical {green) and model (red) error transfer function for the test bench. 
The empirical measurements were taken with an LQG controller which did not model any 
frozen-flow layers. The model was computed using Equation 0 for a 0-layer atmosphere 
(i.e. correcting for static errors only). The transfer function shows the correct roll-off slope 
at low frequency. The model also correctly predicts the controller bandwidth on the real 
system (the lowest frequency where the ETF first crosses 1), though the model underes¬ 
timates the overshoot (the degree to which the system makes frequencies higher than the 
bandwidth worse). 


Figure shows the transfer function measured and computed for an LQG controller which 
does not model any frozen flow turbulence. The LQG controller was implemented in a “DC- 
only” mode (simulating 0 layers of atmospheric turbulence) which has similar performance 
characteristics to an integral controller used in many AO systems. 

The model ETF was computed using Equation Q, but for a single, static layer, and so the 
model behaves like a standard integral controller. The improvements to the test bench allowed 
us to match the roll-off (the slope of the ETE at low frequencies), a feature which is universal 
to all first-order integral controllers, as well as the gain (which sets the left-right position of the 
model in Eigurej^ and bandwidth (the lowest frequency at which the ETF crosses 1.0) between 
the model and the measured transfer function. 

The well matched ETFs show that the system is performing nearly as expected at all spatial 
frequencies, and that more advanced controllers should be stable, with a predictable effect on 
the ETF of the system. 

4.2. Integrating the LQG Controller 

In order to support the LQG controller on our test bench, we made a few modifications to 
the control flow of the test bench. The LQG controller relies on an internal state vector and 
produces an estimate of the entire, open-loop phase at the next timestep. However, both the test 
bench, and the Shane-AO system use the phase error internally and apply the integrators only 
to the commands sent to the deformable mirrors. 

In ShaneAO (and by default on the test bench) the temporal integration is done on the mir¬ 
ror commands directly so that each integrator can consider and provide feedback for actuator 








clipping internally, without affecting the reconstructor or the state of the whole system. This 
is especially important when using multiple deformable mirrors to allow the system to take 
advantage of the high-stroke of the woofer and the many actuators on the tweeter. 

To accommodate this architecture with LQG, we inserted a temporal differentiator into the 
phase reconstructor. The differentiator is applied after the LQG controller is used to predict 
the phase at the next step. Since the commands to each deformable mirror are independently 
integrated, the differentiation at the end of the LQG controller and the integration on each of 
the mirror commands are designed to cancel each other out. This method (differentiating the 
open-loop correction, and applying to the woofer and tweeter separately) is used in GPI 

However, as the individual mirror integrators for the woofer and tweeter DM have differing 
leak settings, and respond differently when actuators clip, the temporally differencing the output 
of the LQG controller is not guaranteed to produce a stable controller. When operating in this 
mode, the effective gain of the system is controlled by the coefficients of the LQG controller, 
while the individual command integrators are set to have a gain of 1, but with their usual leak 
terms to prevent actuator saturation and windup. 

4.3. Demonstrating Stable LQG Control on the MCAO Test Bench 

Once the LQG controller was integrated into the test bench, we were able to demonstrate simple 
stable LQG control. The LQG controller achieved comparable RMS residual wavefront error 
(35.0 nm RMS residual wavefront error from 180.0 nm RMS total wavefront error from the 
atmosphere) to the original integral controller over the course of 4,000 iterations. The far field 
Strehl ratio remained comparably stable during this period. The controller is stable for a variable 
frame delay of 1, 2, or 3 frames (higher delays were not tested). 

Stability was easily maintained even when a wind suppression model was applied with no 
wind translation in the system. The decrease in performance due to an incorrect frozen flow 
assumption was minimal, causing an RMS residual wavefront error of 50 nm. 


4.4. Applying Fourier Wind Identification on the MCAO Test Bench 

We measured the Error Transfer Functions of the test bench with a fast artificial wind speed 
(0.25 subapertures/timestep). The system was run for 2048 iterations in both closed loop, and 
loop open (no correction) mode. Using the closed loop telemetry data, we were able to see 
evidence of frozen flow turbulence at every measured spatial frequency. Figure shows the 
temporal PSD for three Fourier modes (k = 3,6,9 I = 0); each shows clear evidence for frozen 
flow turbulence. 

Given the indication of frozen flow turbulence in the PSDs from the test bench, we were 
able to construct a wind likelihood metric from test bench data showing that the method de¬ 
scribed in Section o can successfully retrieve the input frozen flow wind velocity from the 
system telemetry. Figure shows the identification of the 0.25 subapertures/timestep wind in 
test bench telemetry. 

The identified wind vector (or a-priori known wind vector(s)) can be provided to the Alter 
generation algorithm described in Section 3.2 and derived in m and El which will place a 
notch in the error transfer function for each Fourier mode to remove the temporal peak due to 
the frozen flow layer. The ETF for k = 6, 1 = 0 will be notched near —0.042 timesteps^* and 
will reduce the large power peak visible in Figurej^for that mode. 


4.5. Applying Predictive Fourier Control to the MCAO Test Bench 

To test Predictive Fourier Control, we generated a state-space model for a known windspeed, 
assuming that the wind impacted all Fourier modes. The Alter was generated for wind speeds 
of 0.05, 0.125, and 0.25 subapertures/timestep, and assumed a two frame delay to match our 



system design (see Section]^. To generate the filter, we assumed perfect frozen flow wind was 
detected in all Fourier modes. This assumes that modes where the Fourier Wind Identification 
scheme (outlined in Section [3T| i did not detect peaks due to frozen flow turbulence are never¬ 
theless worth including in the filter. For the test bench, this proved to be the correct assumption, 
but it should be tested on-sky as well. 

We compared the Predictive Fourier Control Kalman filter to a standard integral controller 
in the Fourier Domain, with a gain of 0.25 and a leak of 0.99. The RMS residual wavefront 
error for both the Kalman Filter and the standard Integral Controller is shown in Figure 
The RMS total wavefront error (with no correction) is 180.0 nm. The integral controller has a 
35.0 nm RMS residual wavefront error. Predictive Fourier Control improves the RMS residual 
wavefront error by a factor of 3 to 11.2 nm. 
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Fig. 6. The RMS residual wavefront error for both the Predictive Fourier Controller 
Kalman filter {blue) and a standard integral controller {green) over 4096 iterations. For 
the standard integral, the RMS residual wavefront error was 35.0 nm, and for the Predictive 
Fourier Controller, it was 11.2 nm. The large spikes which occur every 640 iterations occur 
when the phase plate reaches the end of its movement range and must be reset to its initial 
position. With no correction, the RMS total wavefront error is 180.0 nm. 


The measured Error Transfer Function of the test bench when using Predictive Fourier Con¬ 
trol is shown in Figure]^ The ETF was measured using the technique described in Section 
The red dotted line in Eigure[7]shows the expected temporal frequency of the wind for Eourier 
mode k = 6, I = 0 from Equation ([^, ft = —0.042 timesteps^^ A clear “notch” in the PEC 
Kalman filter transfer function is visible at this frequency, where the error transfer function is 
reduced from 0.5 to 10^^. This behavior serves to remove the wind peaks visible in Eigurej^ 
and so eliminates the dominant effect of frozen flow turbulence. 

The elimination of these wind peaks can be seen in the periodogram generated from the 
PEC Kalman filter. The power spectra for three Eourier modes {k = 3,6,9 I = 0) are shown 
in Eigure which lacks the signature peaks from frozen flow that are visible in Eigure 
The figure shows the power spectrum with the active closed loop. Along with the expected 
suppression of power around 0 timesteps^*, there is clear evidence that the Predictive Eourier 
Control Kalman filter suppressed power at —0.042 timesteps^^ Additionally, correction can be 
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Fig. 7. An empirical {blue) and model {red) error transfer function for the test bench 
using Predictive Fourier Control (PFC). For reference, we have plotted a measured trans¬ 
fer function using a standard integral controller {green). The wind notch can be seen at 
—0.042 timesteps^ *. 


applied at or beyond the controller bandwidth, as can be seen by the correction applied tok = 9, 
1 = 0 in Figure § 

4.6. Far Field Performance Improvement 

To correctly predict the potential performance gains of Predictive Fourier Control on a working 
AO system, we also measured far-held images of our light source, which produced an aber¬ 
rated and corrected point spread function (PSF). The PSF was measured every 50 frames (due 
to the slower readout speed of the PSF camera) and then stacked to hnd the effective far-held 
image during closed-loop operation. Using the integral controller, the test bench produced a 
Strehl of 0.479. With the predictive controller, the test bench produced a Strehl of 0.520. These 
measurements contain far-held images from multiple pupil crossings during each test run, and 
were reliably reproduced over multiple runs. Figure shows the radial prohle of the PSF for 
each case. Predictive Fourier Control improves the far-held Strehl performance by approxi¬ 
mately 10%. Gains on higher-Strehl systems should continue to be around 10%. 


4. 7. Comparison of Error Terms 


We used the error budget described in Section 2.2 to compare the performance of the test 
bench when using a standard integral controller, and our Predictive Fourier Control algorithm. 
A complete comparison of the error budgets for each controller is presented in Table 

Using only the standard integral controller, we conhrmed that the test bench can reproduce 
the expected RMS residual wavefront error and expected far-held performance. The test bench 
over performed slightly (by 7 nm) in RMS residual wavefront error, and underperformed in 
the far-held, which we took to be caused by additional non-common-path errors not captured 
by our image-sharpening routine, which, by inference from the performance of the standard 
integral controller, must be 44.0 nm. 
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Fig. 8. Same as Figure]^ but with the Predictive Fourier Control Kalman filter turned on. 
The frozen-flow peaks visible in Figure |^have been suppressed across all Fourier modes. 
The non-zero peak positions are set by Equation The peaks at very high frequencies 
are spatial aliases of the frozen-flow wind. 



Pixets from PSF Center 


Fig. 9. Radial profile of the Point Spread Function of the test bench showing the perfor¬ 
mance improvements between a standard integral controller (green) and Predictive Fourier 
Control (blue). The Strehl improved from 0.479 to 0.520, a relative 10% improvement in 
performance under otherwise identical conditions. 






















Then, comparing the performance of the test bench when using Predictive Fourier Control, 
we found the RMS residual wavefront error to be reduced by 80% of the total time delay error. 
That the Predictive Fourier Controller did not remove 100% of the time delay error is expected, 
as the Predictive Fourier Control transfer function cannot perfectly match the shape of the wind 
peaks in our simple model. This 80% suppression is consistent between the residual wavefront 
measurements and the far-held Strehl measurements, reinforcing that Predictive Fourier Control 
should provide signihcant gains for on-sky adaptive optics systems. 


Table 2. Comparison of Error Budget Terms between Predictive and Standard controllers. 
Terms that are not applicable to the wavefront sensor have been removed. From the final 
performance of the system, we estimate that about 20% of the time delay error remains 
when using PEC. 


Error Term 

Experiment and Controller 

Measurement Source 

WES 

Ear-Eield 

WES 

Ear-Eield 

Reconstructor 

Standard 

Standard 

PEC 

PEC 

Wavefront Sensing 

10.0 nm 

10.0 nm 

10.00 nm 

10.00 nm 

Mirror Fitting 

- 

8.8 nm 

- 

8.8 nm 

Spatial Aliasing 

- 

3.5 nm 

- 

3.5 nm 

Time Delay 

40.5 nm 

40.5 nm 

8.1 nm 

8.1 nm 

Static Uncorrectable 

— 

49.5 nm 

— 

49.5 nm 

Calibration 

- 

44.0 nm 

- 

44.0 nm 

RMS Wavefront Error 

41.7 nm 

78.8 nm 

12.9 nm 

68.1 nm 

Predicted Strehl Ratio 

- 

0.471 

- 

0.522 

Measured RMS Wavefront Error 

35.0 nm 

77.1 nm 

11.2 nm 

68.5 nm 

Measured Strehl Ratio 

- 

0.479 

- 

0.520 


4.8. Performance with Mismatched Filters 

To check the stability of the Predictive Fourier Controller, we ran the test bench using filters 
generated for known incorrect wind speeds. The filters were generated in the same manner as in 
Section [43] but with the incorrect wind direction or magnitude for each test. Using a filter with 
an incorrect magnitude, or an incorrect direction up to 45° from the true wind direction resulted 
in stable correction, with performance degrading to that of the standard integral controller. Wind 
vectors that were significantly far from the true wind on the test bench caused instabilities to 
appear in specific Fourier modes. 

That the filters are stable for minor offsets of the wind vector is unsurprising. Filters with 
I a I < 1 should always be stable El. This suggests that the instabilities we observed with large 
wind-vector offsets are unstable not due to the controller alone, but due to the interaction of the 
controller and the test bench. Future experiments will investigate the specific modal instabilities 
found for large deviations in the wind vector. 

5. Discussion and Future Directions 

This paper demonstrates the feasibility of “Predictive Fourier Control”, but leaves a true char¬ 
acterization of system performance improvements and stability to future work. 

The LQG filter was able to correct for frozen flow turbulence using our Fourier mode state- 
space model. Correction was applied in some modes beyond the controller bandwidth (see 
Figure[7]i, and improved performance at very high temporal frequencies. 








Although we have performed cursory checks of the LQG filter stability using the test bench, 
several questions remain. 

We plan to explore the stability of the LQG filter during real-world operation, where wind 
vectors can vary during the application of the filter. Although the filter appears stable for mod¬ 
erate changes in the wind vector, we will require an on-sky application to test the performance 
gains and losses that occur due to a naturally varying wind vector. 

As the LQG filter maintains a state-space model of the system, it has some “burn-in” time 
during which it must build up a “memory” of the wavefront error. We will characterize, on-sky, 
the effect of this “bum-in” time during filter transitions, and use this to optimize the trade-off 
between the speed of the Fourier Wind Identification loop (the bottom loop in Figure]^ and the 
penalty incurred by changing the LQG filter in real time. 

The PFC algorithm can also be extended to eliminate wind aliasing features, which appear 
as similar planes in Fourier space to those described by frozen flow. In the same way that we 
can reduce delay error in the presence of frozen flow turbulence, we can use the LQG controller 
with the PFC algorithm to correct for the effects of wind aliasing when it is detected in system 
telemetry ll20l . 

Finally, once we have demonstrated the PFC algorithm on-sky with ShaneAO, we intend to 
adapt the LQG controller to account for actuator clipping and the woofer-tweeter architecture. A 
full controller would penalize actuators which are close to clipping within the Kalman filtering 
framework, and so would provide an end-to-end optimal control that would be less reliant on 
the absolute calibration of the individual woofer and tweeter mirrors. 

6. Conclusions 

This paper presents experimental verification of a computationally efficient method for remov¬ 
ing the effects of frozen-flow turbulence from the residual errors in an adaptive optics system. 
Predictive Fourier Control has the potential to reduce the effect of time delay errors, and to 
correct for specifically modeled and known errors beyond the bandwidth of the adaptive optics 
system. 

We demonstrate that in the presence of wind-blown frozen flow turbulence Predictive Fourier 
Control is a feasible method for gaining 10% improvement in an AO system in the far-fleld per¬ 
formance and a factor of 3 in the RMS residual wavefront error. The laboratory experiments 
suggest that we will see significant gains in AO system performance on-sky. This level of perfor¬ 
mance improvement is promising for an algorithm which only requires a software adjustment, 
and which will not add significantly to the computational complexity of an existing adaptive 
optics system. 
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